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1.  INTRODUCTION 


Future  radar  signal  processing  techniques  must  use  doppler  and  polarization  processing 
to  discriminate  between  targets  and  interfering  clutter  and  chaff.  Target  and  clutter 
scatterers  have  characteristic  range,  doppler  and  polarization  properties  which  provide  for 
enhanced  target  detection  performance  compared  to  conventional  single  channel  matched 
filter  receivers.  This  is  due  to  the  fact  that  additional  target  and  background  discriminants 
are  available  in  dual  channel  polarization  sensitive  systems. 

The  scattering  centers  of  radar  targets  and  clutter  often  have  preferred  spatial 
orientations  and  thus  selectively  scatter  radar  signals  according  to  the  orientation  of  the 
incident  electric  field.  In  order  to  accurately  model  radar  targets  or  clutter,  it  is  necessary 
to  consider  their  polarization  properties  in  terms  of  a  complex  scattering  matrix,  which 
represents  a  scatterer's  response  to  two  orthogonal  incident  polarizations.  By  considering 
time  and  range  varying  matrix  elements,  the  scattering  matrix  provides  a  means  for 
jointly  modeling  a  scatterer's  doppler  and  range  dependent  polarization  characteristics. 

This  paper  provides  a  brief  review  of  basic  radar  polarization  theory  [1],  including  a 
discussion  of  the  use  of  the  polarization  scattering  matrix  in  systems  analysis.  Theoretical 
stochastic  doppler  and  polarization  target  and  clutter  models  are  also  discussed,  and  the 
optimum  receiver  is  formulated  [2].  Two  forms  of  the  Optimum  Matrix  Method  [3] 
approximation  to  the  optimum  receiver  are  presented  and  evaluated  using  experimental 
and  simulated  results,  illustrating  receiver  performance  for  the  detection  of  desired  targets 
in  clutter. 

Polarization  discrimination  can  be  exploited  via  three  significantly  different  concepts 
which  may  be  introduced  by  analogy  to  doppler  processing.  In  the  doppler  domain  the 
three  concepts  sire  embodied  as  delay  line  cancellers,  optimum  doppler  filter  banks  with 
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prespecified  transmit  waveforms,  and  doppler  filter  banks  with  adaptive  transmit 
waveforms. 

Thu3,  for  a  delay  line  canceller,  the  design  procedure  is  to  build  a  receive  filter  that 
has  a  null  at  the  center  frequency  of  the  clutter  spectrum.  A  priori  knowledge  of  the  target 
doppler  characteristics  are  not  factored  into  the  design  of  the  delay  line  canceller. 
However,  the  design  of  an  optimum  doppler  filter  bank  depends  upon  both  the  target's 
doppler  frequency  and  the  clutter  spectrum.  Specifically,  each  receive  filter  peaks  at  a 
different  preassigned  doppler  frequency  corresponding  to  an  expected  target  doppler. 
Although  each  filter  of  the  bank  has  high  response  only  to  a  small  range  of  target  doppler, 
the  entire  bank  has  high  response  to  all  possible  target  doppler  frequencies.  The  preceding 
discussion  of  the  two  processors  has  assumed  that  the  transmit  waveform  consists  of  a 
constant  interpulse  period  pulse  train  where  each  pulse  has  the  same  amplitude  and  phase. 
Thus,  as  an  example,  a  two  pulse  canceller  subtracts  and  cancels  successive  returns  from  a 
zero  velocity  point  target  since  the  returns  have  equal  amplitude  and  phase.  However, 
they  have  the  same  amplitude  and  phase  only  when  the  two  sequential  transmitted  pulses 
axe  identical.  The  third  type  of  doppler  processor  adaptively  controls  the  transmitted 
signal  characteristics  to  further  improve  performance.  When  the  amplitude  and  phase  of 
each  transmitted  pulse  can  be  varied,  the  receiver  may  be  designed  to  account  for  this 
particular  waveform,  and  the  resulting  target  doppler  and  clutter  spectrum.  It  has  been 
shown  that  superior  target  detection  performance  is  obtained  via  joint  optimization  of  the 
transmit  waveform  together  with  the  receive  filter,  as  compared  to  optimization  of  the 
receive  filter  alone.  As  the  joint  optimization  depends  upon  the  clutter/chaff  doppler 
characteristics,  the  technique  is  doubly  adaptive.  Both  the  transmit  waveform  and  the 
receive  filter  are  adaptive  to  optimize  the  doppler  discrimination  between  target  and 
clutter. 
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Similarly  for  the  polarization  domain,  there  are  canceller  type  processors  that  are 
designed  to  have  zero  response  at  the  average  polarization  of  the  clutter /chaff  and  whose 
design  does  not  incorporate  any  a  prion  knowledge  of  the  polarization  characteristics  of  the 
targets.  In  addition,  there  is  an  adaptive  receiver  to  enhance  the  discrimination  between 
target  and  clutter/chaff  when  the  transmit  polarization  is  prespecified.  These  two 
techniques  are  compatible  with  doppler  processing.  Thus,  polarization  processing  and 
adaptive  or  non-adaptive  doppler  processing  may  be  cascaded  to  gain  more  effective 
discrimination  between  targets  and  clutter/chaff.  Finally,  aui  adaptive  technique  can  be 
configured  with  both  the  transmit  polarization  and  receive  polarization  adaptable.  This  can 
be  configured  as  dual  domain  adaptivity  by  simultaneously  adapting  in  the  polarization 
and  doppler  domains .  This  last  procedure  is  the  most  powerful  way  of  using  the  two 
domains  for  discriminating  between  targets  and  clutter  /chaff. 

The  goal  is  to  improve  target  detectability  and  to  determine  the  benefits  that  can  be 
accrued  when  the  polarization  and  doppler  domains  are  fully  exploited.  In  order  to 
determine  this,  it  was  necessary  to  assemble  data  to  establish  the  characteristics  of  targets 
and  clutter/chaff.  This  data  base  was  used  to  evaluate  the  absolute  and  relative  target 
detection  capabilities  of  a  number  of  signal  processing  algorithms.  A  particularly  successful 
algorithm  was  a  matrix  method  of  jointly  optimizing  the  transmit  waveform  polarization, 
the  receive  polarization,  and  the  doppler  filter.  This  technique  outperformed  all  other 
approaches  studied. 

The  matrix  approach  is  a  computationally  complex  procedure  which  involves 
simultaneous  optimization  of  the  joint  polarization/doppler  domain.  A  second,  simplified 
technique,  is  one  that  uses  sequential  optimization  of  the  disjoint  polarization  and  doppler 
domains.  Thus,  polarization  optimization  can  be  accomplished  by  using  a  canceller  concept 
based  on  either  polarization  nulling,  polarization  eigenanalysis,  or  receive  polarization 
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orthogonality.  After  adaptive  polarization  processing,  the  doppler  filtering  may  be 
performed  subsequently.  Other  processors  of  interest  for  performance  comparison  are  baaed 
upon  fixed  polarization  systems  using  non-coherent  or  coherent  doppler  processing.  These 
provide  a  reference  to  determine  the  additional  processing  gam  obtained  by  the  optimun 
combined  polarization— doppler  processor  over  conventional  radar  processing  methods. 

The  radar  backscatter  environment  has  stochastic  characteristics  associated  with 
doppler,  range,  and  electromagnetic  scattering  phenomena.  The  environment  consists  of 
target  and  clutter  scatterers,  each  having  characteristic  doppler,  range  and  polarization 
properties.  In  optimization  theory,  Van  Trees  [1]  uses  analytical  approaches  to  solve  a  set 
of  integral  equation.'  for  the  design  of  a  radar  system  optimized  over  a  set  of  radar 
environmental  factors.  This  factors  include  the  effects  of  the  clutter  scattering  function  on 
the  incident  transmitted  waveform.  Such  a  scattering  function  is  directly  related  to  an 
autocovariance  function  of  the  noise  resulting  from  the  backscattered  clutter  with  the 
transmitted  incident  waveform  factored  out.  Such  a  function  for  a  single-channel  system 
incorporates  all  the  statistics  one  needs  to  know  about  the  doppler  and  range  spread 
characteristics  of  the  environment  in  order  to  optimize  the  radar.  For  a  multiple  channel 
system,  the  transmit  waveform  a  a  vector,  and  the  scattering  function  is  a  tensor.  This 
has  been  indicated  by  several  investigators  [2,3]  where  the  parameters  within  the  tensor 
were  time  averaged,  but  doppler  and  range  spread  correlation  were  not  considered.  In  the 
latter  [3],  the  author  indicates  that  the  benefits  can  be  obtained  with  respect  to  the  target 
discrimination  by  controlling  the  transmitted  polarization  waveform  and  receiver  in  such  a 
manner  as  to  take  advantage  of  the  degree  of  coherence  of  the  signals  in  the  receiver 
channels.  There  have  been  many  papers  on  these  subjects.  Stochastic  models  [4] 
describing  the  spectral  and  polarization  characteristics  of  radar  echoes  from  chaff  clouds 
consisting  of  rotating  dipoles  with  completely  random  or  preferred  orientations  have  been 
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used.  Other  papers  [5,6]  include  modeling  and  polarization  processing  with  optimum 
transmit  waveform  vectors  for  the  enhancement  of  target  detection  in  clutter.  These 
approaches  utilize  radar  scattering  matrix  properties  in  the  design  of  a  two-channel  vector 
receiver  which  discriminates  between  the  scattering  matrix  range— doppler  spread 
characteristics  of  the  target  and  those  of  the  clutter.  Computer  simulations  were 
performed  [5]  to  verify  some  simplified  mathematical  models. 

2  THEORY 

This  section  provides  the  analytical  development,  basic  mathematical  expressions, 
and  their  justification,  derived  from  stochastic  processes  and  optimum  filter  theory.  The 
goal  13  the  design  of  the  optimum  system,  making  use  of  the  two  channel  polarization  and 
range-doppler  spread  scatterer  properties.  This  aproach  is  an  extension  of  the  development 
presented  by  Van  Trees  [1]  for  the  single  channel  case. 

A  radar  backscatter  environment  can  be  represented  by  a  statistically  time  varying 
and  range  varying  polarization  scattering  matrix  [7].  The  model  is  general  and  can  be  used 
for  the  clutter  as  well  as  the  target.  Once  the  targets  are  specified,  the  radar  design  can  be 
optimized.  Performance  of  the  system  can  also  be  evaluated.  The  procedures  for  system 
design  and  performance  evaluation  axe  extensions  of  widely  accepted  optimum  filter  design 
theory.  The  primary  difference  is  that  the  target  model  to  which  the  procedures  axe  applied 
is  a  scattering  matrix  random  process.  In  addition  to  serving  as  design  criteria  for  the 
radar  system,  such  a  model  also  designates  the  required  field  measurements  for  completely 
specifying  such  a  scatterer.  As  a  consequence  of  this  model,  an  approach  for  the 
generalized  ambiguity  function  of  multiple  channel  systems  is  introduced.  These 
discriminants  between  the  target  and  the  clutter  will  provide  a  new  dimension  for 
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improving  target  detection  in  a  real-world  clutter  environment. 

First,  we  review  target  modeling  and  system  optimization  for  the  single  channel  radar 
system.  Using  the  review  as  a  basis,  we  will  extend  the  treatment  of  the  multiple  channel 
system  specialized  for  the  polarization  case.  The  procedures  for  deriving  the  target  model, 
the  optimum  transmit  waveform,  and  the  optimum  receiver  for  the  single  channel  case,  or 
single  polarization  sense,  axe  given  as  follows: 

The  signal  received  by  the  radar,  r(t),  from  a  slowly  fluctuating  point  target  in  a  clutter 
environment  is 

r(t)  =  b  fd(t)  +  nc(t)  +  w(t)  ,  Hj  (1) 

r(t)  =  nc(t)  +  w(t)  ,  HQ  (2) 

where: 

hypothesis  target  present,  Hq  target  absent 
f^(t)  is  the  complex  range  delayed  doppler  shifted  replica 
of  the  transmit  waveform  f(t) 

b  is  a  complex  random  variable  and  represents  the  target 

backscatter,  propagation  losses,  and  antenna  responses 
nc(t)  is  the  received  signal  from  the  clutter  and  is  a  complex 
random  process 

w(t)  is  the  additive  white  noise  in  the  receiver. 

Tb*  tilde  ~  designates  the  complex  low-pass  equivalent  of  the  Tadar  signal. 

The  received  signal  from  the  clutter  is  the  convolution  of  the  clutter  with  the  transmit 
waveform.  Convolution  is  with  respect  to  the  range  variable  X. 

n.(t)  =  /ij  f(t-A)b(t-£,  A)<U 

-CD 
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Et  is  the  energy  in  the  transmitted  waveform. 

n  (t)  is  a  zero-mean  complex  Gaussian  random  process  with  the  covariance  function: 


n 


/ 


CD 


K  ~  (t,u)  =  Ef  J  f(t  -  A)  KDR{t  -  u}  f  (u  -  A)  dA 


=  Jj  f(t  -  u)  SpR{i^,A}  f  (u-A)  e^  u)dt/  dA 
—00 


(4) 


(5) 


where  the  correlation  function,  Kj-)R{t,A},  i®  a  two-variable  function  that  depends  on  the 
reflective  properties  of  the  target,  and  Sj-jR{^,A}  represents  the  spectrum  of  the  process 
b(t,A)  and  is  called  the  scattering  function  of  the  clutter.  The  two  functions  are  a  Fourier 
Transform  pair.  The  white  noise  w(t)  is  likewise  a  complex  Gaussian  random  process  but 
with  covariance  function  K  *  (t,u)  =  Nq  £(t-u). 

The  conventional  receiver  is  designed  for  detecting  a  target  in  white  noise  only.  We 
compute  the  sufficient  statistic  £  : 


r(t)  f  (t  -  t^)  e-^  wd  ^  dt 

-00 

and  make  the  threshold  comparison  test: 

H1 

I  l  I2  |  7 


(6) 


H0 

The  output  (S/N)  is  given  by: 


(7) 


u. 


*  /  No 

1  +  'r 


E  /  N 
r  '  o 


1  +  (Et/No)  jj  *  dv  dA  Sdr{*/,A}  0{rd- A,  ^fd} 


(8) 


(9) 
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/>r  lepresents  the  degradation  due  to  the  clutter 
f(t)  is  the  transmit  waveform  which  is  used  here  as  a 
modulation  or  filter  function 
Tj  is  the  target  round  trip  delay 
Wj  is  the  target  doppler  frequency 

7  is  a  predetermined  threshold 

1?.  is  the  average  energy  received  from  the  target 

is  the  transmit  energy 
Nq  is  the  white  noise  energy  (kT) 

0{A,i/}  is  the  waveform  ambiguity  function 


The  optimum  single-channel  receiver  is  designed  for  detecting  a  target  in  clutter  as 
well  as  white  noise.  The  optimum  receiver  computes 

r  -  f  ’  r*(t)  g'«  *  (10) 

—05 

where  g(  • )  satisfies  the  integral  equation 

f (t  -  Tj)  e  d  jfjf  f(t  -  A)  Kpp(t-ut  A)  f  (u-A)  g(u)  du  dA 

“CO 

+  N0  g(t)  (11) 

The  function  g( •)  is  a  modulation  or  filter  fimction  and  is  used  in  place  of  f( •)  when 


colored  noise  is  present.  The  performance  is  obtained  by  evaluating 
T 

A  =  T)  Jjf  ^  f(t)  Q-  (t,u)  f  (u)  dt  du 

r  *'*'  T  n 

i 

T 

A  =  E  f(t)  i*(t)  dt 

where  the  inverse  Kernel  Q-(t,u)  satisfies 

n 


(12) 

(13) 
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/ 


V- 


K-  (t,.)  +  N 


Q~(z,u)  dz  =  $(t-u) 
n 


(14) 


Then  the  false  alarm  and  detection  probability  are  related  by 


PFA  =  (pD>1+4  <I5> 

For  certain  special  conditions  of  the  scattering  function  we  cam  simplify  equations  (10) 

and  (11)  further.  The  optimal  waveform,  fQ(t),  must  satisfy  [1]  the  following  integral 

equation 


T 

/  f  hou(t,u)  ~fo^  du  +  AE  *  ^  ~  XB  V)  =  0  (16) 

i 

where  hQu(t,u)  is  the  optimum  unrealizable  filter  satisfying  the  equation 

N^fcu)  +  (t,x)dx  =  K~  (t,u) 

T-  n  n 

i  c  c 

T;  <  t,u  <  Tf  (17) 

Ag  and  Ag  are  LaGrange  multipliers  with  an  energy  and  bandwidth  constraint. 

When  we  transmit  and  receive  over  two  channels  as  we  do  in  the  dual  polarization 
case,  the  received  signals  given  in  Equations  (1)  and  (2)  become  two-element  vectors 
r(t)  =  bfd(t)  +  nc(t)  +  w(t)  Hx  (18) 

r(t)  =  n.  w(t)  Hq  (19) 

where  boldface  variables  denote  vectors  of  the  polarization  components.  The  vector 
received  signal,  r(t),  contains  a  component  from  each  receive  polarization  channel. 

The  scattering  matrix,  b(t),  contains  four  components:  the  vertically  polarized  scattering  of 
the  incident  vertically  polarized  energy  (VV),  the  cross  polarized  components  (VH  and 
HV),  and  the  horizontal  scattering  component  (HH).  The  scattering  is  modeled  as  a  slowly 
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fluctuating  point  target  with  zero  mean  complex  Gaussian  random  elements.  In  general, 
the  scattering  components  byg  and  bgy  may  be  assumed  to  be  equal.  f^(t)  is  the  time 
delayed  doppler  shifted  replica  of  the  transmit  waveform  vector  with  transmit  energy  E^. 
The  receiver  noise,  w(t),  is  a  vector  with  white  noise  components,  Wy(t)  and  Wg(t),  which 
are  independent  for  the  two  polarization  channels.  The  colored  noise  vector  due  to 
backscatter  from  clutter,  n  (t),  also  contains  dual  polarization  components  which  are 
obtained  by  convolving  the  range-doppler  variant  scattering  matrix  of  the  clutter  process, 


b(t-|,A)4 


"by(t  -  A) 

A 

'byy(t-'J)  A)  byg(t_'J1^) 

>HV^  ■  7’  A)  ~  ■?’ 

(20) 


with  the  transmit  vector, 


f(t)  = 


fyOO 

J’hW- 


(21) 


which  is  normalized  to  unit  energy.  The  scattering  matrix  b(t— A),  a  zero  mean  Gaussian 
random  process,  interacts  with  the  transmit  waveform  resulting  in 


nc(t)  = 


v/*) 

ncHW 


=  Je. 


b(t  -  4*  A)*M)  dA 

—00 


(22) 


This  expression  represents  the  colored  noise  vector  due  to  baclcscatter  from  the  clutter, 
where  t  is  time  and  A  is  radar  range  expressed  in  units  of  time.  It  can  be  shown  that  the 
clutter  vector  covariance  function  is  a  matrix 
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statistically  independent  and  that  the  return  from  each  interval  ia  a  sample  vector  function 
of  a  stationary  zero-mear  Gaussian  random  process, 


K;  (t,u)  =  Et  f  ”  fT(t  -  A)  KDR{t-u,  A}  ?(u  -  A)  dA 


--  ®  j2n/(t-u) 

fi(t_  A)  SDR^’A^f  (u_  A)  e  dvdX 

—oo 


(24) 

(25) 


The  tensor  correlation  function  KjjR{r,A}  is  a  two-variable  fourth-rank  tensor  that 
depends  on  the  reflective  properties  of  the  clutter.  The  ‘DR1  subscript 
denotes  that  the  clutter  is  doubly  spread  in  both  doppler  and  range.  The  function 
,A}  is  a  two-variable,  fourth-rank  tensor  representing  the  spectrum  of  the  process, 
and  is  related  to  Kj-jR(t,A)  by  the  Fourier  transform 


00  - 


_  —  flm/T 

^DR^7’^  =  f  e  dv 


(26) 


It  can  be  called  the  tensor  scattering  function  of  the  process  b(r,A).  KpR(r,A) 
provides  16  different  elements  (discriminants)  when  one  considers  the  statistical  behavior 
of  the  polarization  random  process  scattering  matrix.  Assuming  the  scattering  matrix  is  a 
non-negative  Hermitian  process,  then  the  number  of  independent  elements  reduces  to  10, 
which  still  provides  a  wide  selection  of  polarization/doppler  discriminants.  These  tensor 
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element  functions  completely  describe  the  clutter  irrespective  of  the  transmit  waveform 
and  receiver  design.  They  can  be  used  in  the  waveform  and  receiver  optimization  as  well 
as  the  performance  equations  in  the  same  sense  that  the  scalar  correlation  and  scattering 
functions  sure  used  in  the  design  equations  of  the  single  channel  system. 

The  conventional  receiver  is  one  which  is  optimized  for  a  signal  corrupted  by  white 
noise  instead  of  the  clutter  colored  noise.  In  this  case,  the  optimum  filter  is  a  vector 
matched  filter.  The  sufficient  statistic  computed  by  the  conventional  receiver  is, 

e  -  /  “  iT(t)  f*(t  -  Td)  e  J"d*  dt  (27) 

-co 

and  is  compared  to  a  threshold  as  in  equation  (7) .  The  performance  degradation  is  given 
by  equation  (9).  For  the  vector  case,  we  leave  the  ambiguity  function  in  tensor  /integral 
form. 


/>r(Td>u'd)  ~  TT  f  f  f  f  ^t_rd)  *I'(t“A)  e 


j(w-wd)t  - 

u- 


DR 


(M) 


~  -j(aH-a/j)u 

f  (u-A)  f(u-rd)  e  dt  du  dA  di/  (28) 


Equation  (28)  introduces  another  model,  that  of  the  ambiguity  function  for  the 
multi-channel  case.  It  takes  on  a  form  suggested  by  the  integral  of  the  outer  product  of  all 
the  transmit  vectors  f  in  equation  (28) . 

In  equation  (9) ,  the  integral  contained  the  waveform  ambiguity  function 
8  {t^-A,  t'-t'j}.  Equation  (28)  generalizes  this  function  for  the  multi-channel  case  and 
contains  the  necessary  functions  and  integrads  to  constitute  an  ambiguity  function. 
However,  it  involves  outer  products  with  waveform  vectors.  The  optimum  receiver 
computes 
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(29) 


a  / "  ;T(t)  r  (t)  ^ 

— tu 

for  the  vector  case,  and  makes  a  threshold  comparison  as  cited  in  equation  (7). 
A  block  diagram  of  the  optimum  receiver  is  given  in  Figure  1. 


V'1 


Figure  1.  Polarization  Modulation  Receiver 
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g(  •)  satisfies  the  matrix  equation: 


fd(t)  =  /  ^  (t,u)  g(u)  du  +  Nq  g(t) 

-CD  C 


=  ff°  fT(t-A)  KDR(t-u,  A)  f  *(u-A)  iu)  du  dA  +  Nog(t) 


The  performance  (output  signal  to  noise  ratio)  of  the  optimum  receiver  is  specified 


A  ~  ^r  ff  f  ^  (t)U)  *  ^  dt  du 

1 

=  J  dt 


The  average  energy  received  from  the  target  is 


Er  =  E  /  ”  ST(t)  S*(t)  dt  •  Et 
— 00 


where  S(t)  =  (t). 

d 

This  reduces  to 

^r  =  2  (*vv  +  ^Hv)  EtV  +  2  ^HH  +  EtH  + 

^VV.VH  +  °HV,Hh)} 


where  ^yy  18  the  variance  of  the  byy  element  in  the  target  scattering  matrix  and 

2 

<Tyv  yg  is  the  covariance  between  the  byy  and  the  byg  elements  of  the  target 
scattering  matrix.  E^y  and  E^g  are  the  transmit  energy  of  the  vertical  channel  and  the 
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horizontal  channel,  respectively,  p ^  is  the  correlation  coefficient  between  the  vertical  and 
horizontal  waveforms. 

The  quantity  A  for  which  expressions  have  been  given  throughout  this  discussion  is 
used  to  relate  the  detection  and  false  alarm  probability 

PF  =  (PD)1+A  (3«) 

The  above  optimization  procedures  involve  integral  and  tensor  equations  which  are 
very  difficult  to  solve  analytically.  In  general  even  the  scalar  case  is  not  trivial.  Numerical 
methods  applicable  to  computer  approximations  of  these  functions  are  valuable  in 
approaching  the  optimum  performance. 

The  extension  of  the  waveform  design  case  to  the  polarization  channels  involves 
vectors,  matrices,  and  tensors  in  equations  (16)  and  (17).  The  optimum  unrealizable  filter, 
hQU,  is  a  2  x  2  matrix.  The  optimum  waveform,  fQ  is  a  vector,  and  the  autocorrelation 
function,  Knc  is  a  matrix.  It  is  important  to  note  that  the  autocovariance  matrix  contains 
factors  of  the  waveform  vectors  as  defined  in  equation  (24).  Consequently,  the  filter  must 
satisfy  the  equation 


T  CD 

"oVM  =  Et/T7  fl(t-A)  KDR{t-x,A} 

Ti  _ 

[I^x-u)  -  hou(x,u))  dA  dx 


and  the  optimal  waveform  vector  must  satisfy  the  integral  equation 


/Tf  Eou(t'u)  'f0<u>  du  +  AE  "foM  -  AB  =  0 


(37) 


(38) 
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expression  for  the  performance  of  the  extended  processor.  In  the  extended  processor  case, 
the  receiver  contains  vector  channels,  each  having  its  own  filter  or  modulating  function 
determined  from  the  time/range  statistical  properties  of  the  transmit  polarization 
waveform  reacting  with  the  clutter/target  model.  We  also  gave  an  expression  which  an 
optimal  transmit  waveform  must  satisfy  for  the  polarization  case. 

Other  approaches  have  been  taken  for  the  optimal  processing  of  polarimetric  received 
signals.  Mohrhotra  [8]  decomposes  the  return  signal  into  a  complete  orthonormal  (C  O  N.) 
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set  of  eigenvalues.  He  develops  all  the  associated  processing  in  a  similar  manner  as 
developed  in  this  paper  with  an  additional  stage  constituting  signal  and  channel  weights  as 
functions  of  the  eigenvalues.  Each  weight  is  applied  to  the  channel  in  which  the  received 
signal  has  been  correlated  with  its  associated  eigenvector  function. 

An  optimum  receiver  which  correlates  the  received  signal  r(t)  with  g(t),  the  vector 
modulation  function,  requires  exact  knowledge  of  the  scattering  statistics  of  the  clutter  for 
optimality.  Since,  in  practice,  only  a  numerical  approximation  of  the  clutter  covariance 
matrix  is  available,  and  this  information  is  obtained  by  a  limited  time  observation  of  the 
clutter  itself,  a  linearization  of  the  optimum  receiver  suffices  to  compute  detection 
probability  versus  signal-to-noise  ratio  for  a  selected  false  alarm  probability.  A  further 
simplification  is  imposed  by  high  power  transmitter  design,  which  favors  the  use  of 
waveforms  with  uniform  amplitude  characteristics. 

The  receiver/processor  implemented  in  the  RADC  Surveillance  Laboratory,  called  the 
Optimum  Matrix  Method  (OMM),  follows  the  foregoing  approach.  Based  on  joint 
diagonalization  of  the  target  and  clutter  covariance  matrices,  OMM  provides  maximum 
separation  of  the  target  and  clutter  vector  components  in  the  processed  output,  prior  to 
hypothesis  testing.  The  method  includes  transmit  vector  waveform  optimization  derived 
from  estimated  target  and  clutter  statistics.  An  efficient  waveform  design  technique  has 
been  developed  [5]  which  selects  a  most  nearly  optimum  candidate  from  an  assembled 
library  of  constrained  transmit  waveforms.  Four  other  receivers  and  processing  techniques, 
exhibiting  various  degrees  of  polarization  and  doppler  adaptivity,  have  been  evaluated  in 
several  experiments  for  comparison  with  the  OMM,  and  are  detailed  in  section  3.4 
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3.0  APPROACH 


This  section  covers  the  processor  evaluation  which  includes  experimental  as  well  as 
simulated  data  utilized  to  predict  the  detection  performance  of  various  algorithms. 
Performance  is  measured  in  terms  of  probability  of  detection  (PD)  versus  signal-to-noise 
ratio  (SNR)  as  a  function  of  clutter-to-noise  ratio  (CNR).  The  effort  was  divided  into 
three  areas:  data  collection,  data  reduction,  and  processor  evaluation . 

3.1  DATA  COLLECTION 

The  experimental  data  sets  were  obtained  using  the  S-Band,  dual  polarized  tracking 
radar  testbed  located  at  Rome  Air  Development  Center  (RADC),  in  Rome,  New  York. 
The  RF  transmitting  and  receiving  equipment,  and  antenna,  Eire  located  on  a  65  foot  tower 
adjacent  to  the  RADC  Surveillance  Laboratory,  as  shown  in  Figure  2.  The  antenna 
features  Ein  elevation  over  azimuth  pedestal  configuration  which  provides  monopulse 
tracking  and  variable  polarization.  The  main  central  horn  is  duEil  polarized  for 
transmitting  and  receiving.  Tv.  o  independent  S-Bard  channels  are  available  on  trEinsmit 
and  receive.  On  receive,  the  two  orthogonally  polarized  channels  sire  synchronously 
recorded  and  processed  independently.  The  system  has  two  transmitters,  one  for  each 
polarization  channel.  This  approach  has  the  eid vantages  of  independent  frequency 
capability,  polarization  versatility,  and  increased  power,  relative  to  a  single  transmitter 
configuration  Transmit  polsinzation  agility  may  be  achieved  at  low  power  levels  by 
controlling  the  signal  inputs  to  the  two  transmitters.  A  wideband  mode  is  available  for  the 
S-Band  waveform  generation  and  processing  systems  which  allows  signal  bandwidths  of  up 
to  320  MHz  to  be  utilized. 
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The  S-Band  radar  system  parameters  are  as  follows: 


System: 


Frequency 

3.1  -3.7  GHz 

Peak  Power 

170  Kilowatts 

Average  Power 

15  Kilowatts 

Pulse  Width 

Up  to  200  usee  (max) 

PRF 

Up  to  5000  pps  (max) 

Dynamic  Range 

+90dB  in  1MHz  bandwidth 

Data  Channels 

Four:  Two  I,  Two  Q 

Cassegrain  Feed: 

Polarization  Five  Horn  Cassegrain 

Sum  and  Error  Horns  -  Linear 

Isolation  30dB  between  vertical  and 

horizontal 

Tracking  Null  30dB  between  difference  and  sum 

channels 
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Antenna: 


Gain 

Dish  Diameter 
Sidelobes 


Focus  to  Diameter 
Azimuth  Travel 
Elevation  Travel 


42dB  for  main  channel 
23dB  for  monopulse  channel 
20  feet 

Elevation  and  Azimuth,  1st  and 
2nd  sidelobes  20dB  below  main. 
3rd  is  25dB  below  main. 

0.4 

720°  (wrap  limited) 

-5°  to  ld5° 


The  data  collection  activity  consisted  of  recording  dual-polarized  returns  from  targets 
and  clutter  using  the  3-Band  system.  Figure  3  is  the  experiment  block  diagram.  The 
S-Band  antenna  provided  independent  outputs  for  the  horizontal  (H)  and  vertical  (V) 
receive  channels.  Both  channels  were  synchronously  detected  by  the  receivers  whose 
baseband  outputs  were  digitized  using  12  bit  analog  to  digital  converters.  The  A/D 
converter  outputs  were  buffered  and  transcribed  to  VAX-compatible  9-track  digital 
magnetic  tapes  for  analysis  by  software  realizations  of  several  processing  algorithms.  By 
operating  on  recorded  data  from  the  radar  data  base,  a  valid  comparison  was  possible  by 
replaying  flight  test  data,  as  opposed  to  conducting  severed  live  tests  for  the  different 
processors.  Figure  4  shows  the  equipment  configuration  for  the  data  collection  activity. 
The  upper  part  of  the  diagram  shows  the  on-line  system  that  recorded  data  onto  the  high 
bit  rate  (HBR)  instrumentation  recorder,  while  the  lower  part  shows  the  off-line  data 
processing  used  to  transfer  sections  of  data  to  6250  bpi  9-track  magnetic  tapes. 
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The  on-line  system  used  to  gather  data  was  the  existing  Signal  Processing  Laboratory 
(SPL)  S-Band  radar  system  with  the  HBR  recorder  interfaced  to  the  receiver  outputs.  The 
S-Band  antenna  was  controlled  by  hardware  range  and  angle  trackers.  The  transmit 
section  consisted  of  the  SPL  waveform  generator  driving  the  S-Band  exciter,  the  vertical 
and  horizontal  transmitters,  and  the  transmit  polarization  controller. 

The  waveform  generator  produced  a  linear  frequency  modulated  (LFM)  waveform  that 
was  applied  to  both  transmit  channels  simultaneously.  The  transmit  polarization  controller 
operated  the  transmitters  sequentially  so  that  interleaved  horizontally  and  vertically 
polarized  pulses  were  produced. 

The  receive  section  included  the  duplexers,  receivers,  low-pass  filters,  and  analog  to 
digital  converters  for  the  vertical  and  horizontal  channels.  The  converter  outputs  were 
formatted  by  the  2:1  multiplexer  for  input  to  the  HBR  recorder.  This  machine  recorded 
A/D  data  at  a  complex  word  rate  of  4  MHz  continuously  in  real  time. 

Radar  parameters  used  for  all  experiments  were  as  follows. 


Operating  frequency 

3.35  GHz 

Peak  output  power 

150KW,  nominal 

Transmit  PRI 

2  milliseconds 

Waveform  type 

Linear  FM 

Waveform  bandwidth 

2  MHz 

Pulse  duration 

32  microseconds 

Pulse  compression  ratio 

64  to  1 

A/D  converter  resolution 

12  bits 

A/D  converter  sample  rate 

2  MHz 
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The  off-line  system  retrieved  data  from  the  HBR  tape  and  transferred  sections  to  hard 
disk  on  an  HP  A900  computer  in  the  SPL,  then  to  VAX-compatible  9-track  tape.  The 
intermediate  transfer  to  disk  storage  was  required  to  keep  pace  with  the  high  playback  data 
rate  of  the  HBR  recorder.  The  resulting  9-track  tapes  could  then  be  processed  with  no 
limitations  on  algorithm  complexity  and  computation  times. 

The  procedure  used  to  perform  the  data  collection  activity  was  as  follows: 

a.  Configure  the  dual-polarized  S-Band  radar  to  transmit  interleaved  horizontally 
and  vertically  polarized  LFM  pulses  while  receiving  both  polarization  channels 
simultaneously. 

b.  Calibrate  the  radar  in  order  to  achieve  a  match  between  the  horizontal  and  vertical 
receive  channels  in  magnitude  and  phase.  Independently  match  the  H  and  V  transmit 
channels  in  magnitude  and  phase. 

c.  Engage  the  acquisition  radar  system  to  locate  targets  and  clutter  areas  of  interest. 
The  azimuth  and  elevation  readings  initialized  and  cued  the  S-Band  antenna  tracker. 

d.  While  tracking,  record  several  minutes  of  the  radar  data  onto  the  HBR  tape 
recorder. 

e.  Transfer  several  thousand  PRI  records  of  receiver  data  from  the  HBR  tape  to 
VAX-compatible  9-track  tape  via  hard  disk  on  the  HP  A900. 
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3.2  DATA  REDUCTION 


The  primary  function  of  the  data  reduction  activity  was  to  interface  between  the  raw 
A/D  converter  outputs  recorded  on  tape  and  the  processor  evaluation  software.  The 
processor  evaluation  was  based  on  exploiting  the  different  polarization  and  doppler 
characteristics  of  the  target  and  clutter  under  analysis.  These  characteristics  were 
contained  in  a  set  of  six  covariance  sequences,  as  outlined  below,  derived  from  a  particular 
section  of  raw  data  by  the  data  reduction  procedure  and  placed  in  a  FORTRAN  data  file. 
Data  reduction  was  performed  on  a  number  of  sections  of  recorded  returns  to  yield  a 
collection  of  covariance  sequence  files.  During  the  processor  evaluation  activity,  these  files 
were  used  as  inputs  to  the  software  that  calculated  PD  curves  for  the  processors  of  interest. 
Three  examples  are  presented  in  Figure  5. 


SINGLE  POLARIZATION 


VARIABLE  POLARIZATION/SINGLE  RCVR 


VARIABLE  POLAPIZATION/MUI.TIPLE  RCVR 


Figure  5.  Three  Polarization  Receivers/Processors 
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A  second  function  of  data  reduction  was  to  provide  additional  information  for  judging 
the  quality  of  the  return  data.  Normally,  a  section  of  contiguous  pulses  was  processed  from 
the  digital  tape  at  one  time.  The  data  reduction  software  produced  compressed-pulse  plots 
of  the  horizontal  and  vertical  return  magnitudes  versus  range  cell  for  the  first  pulse  in  the 
selected  section.  All  pulses  in  the  section  were  then  pulse-compressed.  A  particular  range 
cell  of  interest  was  chosen  for  analysis  and  the  pulse  compressor  output  values  for  this 
range  cell  in  each  return  pulse  were  stored  for  subsequent  processing.  The  data  reduction 
software  printed  these  peak  values  in  a  table  and  plotted  their  magnitudes. 

The  third  function  of  the  data  reduction  activity  was  to  calculate  and  print  out  vector 
autoregressive  model  coefficients  for  the  process  generating  the  return. 

3.3  PROCESSING  ALGORITHMS 

The  process  used  for  calculating  the  covariance  sequences  begem  by  reading  in  and 
formatting  the  digitized  radar  samples  from  9-track  tape.  The  radar  transmit  schedule  was 
chosen  during  data  collection  such  that  alternating  horizontal  and  vertical  pulses  were 
transmitted.  The  receive  system  recorded  two  independent  polarization  channels 
simultaneously  on  each  PRI.  When  a  horizontal  pulse  was  transmitted,  two  receive  signals 
were  captured,  designated  the  horizontal  transmit/horizontal  receive  (HH),  and  the 
horizontal  transmit/vertical  receive  (HV)  returns.  During  the  next  PRI,  a  vertical  pulse 
was  transmitted,  and  the  two  receive  channels  captured  the  vertical  transmit/horizontal 
receive  (VH),  and  the  vertical  transmit/vertical  receive  (VV)  returns.  The  data  formatting 
consisted  of  sorting  these  four  types  of  returns  into  separate  data  streams. 

The  return  pulses  were  pulse-compressed  with  a  time-bandwidth  product  of  64.  This 
was  done  in  software  using  a  two-transform  technique  with  35  dB  Taylor  weighting  for 
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time  si  delobe  suppression. 


The  six  covariance  sequences  generated  were: 


RHH,HH(n>  =  E<HHW  HH  Vn)}.  1  <  i  <  N  (41) 

RVV,VV(n)  =  E<VV«  VV*(i+n)},  1  <  i  <  N  (42) 

RHV,HvW  =  E<HVM  HV*(i+n)},  1  <  i  <  N  (43) 

RHH,VV(n)  =  E{HH(i)  w’(l+n)},  1  <  l  <  N  (44) 

RHH  Hv(n)  =  E{HHM  1  <  >  <  N  (45) 

RVV,HV<n>  =  E^VV(i)  HV*(i+n)},  1  <  i  <  N  (46) 


where  R(n)  is  the  complex  covariance  sequence  with  the  lag  index  n,  E{  • }  is  the  expected 
value  operator,  and  *  is  the  complex  conjugation  operator.  N  is  the  number  of  pulses 
processed.  The  arrays  HH,  VV,  and  HV  contain  the  complex  values  of  the  compressed 
pulses  at  the  particular  range  cell  selected  for  processing. 

The  cross-polarization  data  array  was  formed  by  interleaving  compressed  pulse  values 
at  a  single  range  cell  from  the  HV  and  VH  data  streams.  This  represents  the  array  of  all 
cross-polarized  returns  from  the  chosen  range  cell.  For  a  64-pulse  window,  there  were  32 
HV  returns,  and  32  VH  returns.  The  compressed  pulse  values  from  each  range  cell  of 
interest  were  put  into  the  HV  array  so  that  the  first  element  was  an  HV  sample,  the  second 
element  a  VH  sample,  etc. 

The  HH  and  W  arrays  were  formed  differently  from  the  cross-pol  array  because  each 
has  only  one-half  the  number  of  samples  as  the  interleaved  HV/VH  array.  The  samples 
filled  the  odd  indexed  elements  of  the  HH  array.  The  VV  array  was  filled  so  that  the  even 
indexed  elements  contained  the  sample  data.  Since  there  were  only  one-half  as  many 
observations  of  the  co-polarized  returns  as  that  of  the  cross-polarized  returns,  the  values  in 


27 


the  HH  and  VV  arrays  were  interpolated.  The  result  was  a  set  of  three  arrays,  one 
containing  cross-polarized  observations,  one  containing  HH  observations,  and  one 
containing  the  VV  observations.  The  cross-pol  and  interpolated  HH  and  VV  arrays  each 
contained  values  representing  observations  at  intervals  of  one  PRI.  All  the  arrays 
contained  the  same  number  of  elements,  equal  to  the  number  of  PRIs  in  the  pulse  window 
to  be  processed. 

The  covariance  sequences  were  estimated  from  these  arrays.  This  was  accomplished 
using  a  4096  point  Fast  Fourier  Transform  and  the  results  normalized  to  zero  doppler 
frequency  before  the  estimates  were  made.  The  samples  were  shifted  so  that  the  doppler 
peak  was  moved  to  *he  first  sample.  The  location  was  determined  by  a  magnitude  peak 
search.  Doppler  normalization  was  required  for  the  target  because  the  processor  evaluation 
software  required  target  data  at  zero  doppler  frequency.  Normalization  was  not  performed 
on  the  clutter  data. 


The  covariance  sequences  were  computed  for  all  six  combinations  of  the  HH,  W,  and 
HV  sample  streams.  They  carried  the  doppler  and  polarization  characterization  of  the 
target /cl utter  needed  for  use  by  the  processor  evaluation  algorithms. 

Given  the  experimentally  derived  covariance  functions  for  each  of  the  scattering 
matrix  elements,  the  stochastic  vector  h(t)  was  formed  by  modeling  the  data  via  a  complex 
vector  autoregressive  process.  The  scattering  model  thus  obtained  provided  the  basis  for 
the  Optimum  Matrix  Method  design.  This  allowed  higher  resolution  estimates  to  be 


performed  on  limited  observation  data,  enhancing  performance. 

An  M^  order  autoregressive  process  modeling  the  stochastic  vector  h(t)  is  defined  as 
follows. 


M  ~  ~ 

h(t)  =  £  [  A,  h(t-kT)  +  u(t)  ]  (47) 

k=  1  K 
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where 

=  3x3  matrix  of  complex  coefficients  to  be  determined 
h(t)  =  3x1  complex  vector  of  covariance  sequences 
= 

-  hw^  • 

u(t)  =  3x1  complex,  zero-mean  white  Gaussian  noise  vector 
T  =  Sampling  period  (radar  PRI) 

M  =  Order  of  AR  model 

The  A^  matrices  are  obtained  by  solving  the  Yule-Walker  equations 

M  -  - 

K(mT)  =  £  A,  K[(m-k)  T]  m  =  1,  2,  3,  ...  M  (48) 

k=l  K 

where 

K(mT)  =  3x3  complex  covariance  matrix 

=  E{  h(t)  h  (t-fmT)  }  (49) 


The  M  autoregressive  model  coefficient  matrices  are  determined  by  solving  equation 
(48)  M  times.  The  K  matrices  are  determined  by  the  data  reduction  procedure: 


K(mT)  -  RHH,HV^m^  RHH,VV^m^ 

RHH,Hv(m)  RHV,HV^m^  ^V.HV^ 

D  / _ \  o  / _ \  r>  / _ \ 


These  are  the  tensor  covariance  elements  of  the  scattering  matrix  developed  in  the 
theoretical  discussion  of  the  previous  section  with  the  assumption  that  the  scatterers  are 
isotropic,  that  is,  byg  = 
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3.4  PROCESSOR  EVALUATION 

Five  different  processing  algorithms  were  implemented  in  software  and  are  briefly 
described  here. 

The  Optimum  Matrix  Method  (OMM)  is  shown  in  block  diagram  form  in  Figure  6. 
The  processor  consists  of  a  weighting  matrix  followed  by  a  maximum  likelihood  receiver. 
The  inputs  to  this  processor  are  the  target  and  clutter  samples  of  the  scattering  matrix. 
The  first  section  of  the  processing  program  computes  the  optimum  transmit  waveform  and 
the  receive  weighting/transform  matrix.  The  second  part  of  the  program  consists  of  an 
iteration  loop  wherein  the  target  and  clutter  samples  are  processed  by  an  optimum 
receiver.  This  includes  transmitting  the  optimum  waveform,  operating  on  the  return  with 
the  weighting/transform  matrix,  and  comparing  the  result  with  a  threshold.  The  number  of 
threshold  crossings  is  then  counted  and  used  in  determining  detection  performance.  A  block 
diagram  of  the  vector  weighting  method  is  shown  in  Figure  7.  The  method  consists  of  two 
parts.  The  first  includes  an  iterative  loop  to  calculate  the  optimum  weighting  vector  then  a 
computation  of  the  signal-to— clutter  ratio  resulting  from  the  application  of  this  optimum 
vector.  The  second  part  computes  the  signal-to-clutter  ratio  for  the  simple  matched  filter 
case.  Plots  of  the  outputs  provide  comparison  of  the  optimum  vector  approach  to  the 
conventional  matched  filter  method. 

The  dual  adaptive  processor  used  for  evaluation  of  the  optimum  matrix  method  is 
described  in  reference  [5].  It  is  implemented  in  two  forms.  The  first  is  the  original  form,  m 
which  both  transmit  and  receive  polarizations  sue  adaptive  to  any  orientation.  The  second, 
designated  the  constrained  optimum  matrix  method  (COMM)  processor,  can  adapt  the 
transmit  polarization  only  to  one  of  sixteen  discrete  polarization  orientations. 
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Figure  6.  Optimum  Matrix  Method  Processing  Algorithm 


Figure  7.  Block  Diagram  of  Vector  Weighting  Simulation 
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The  next  technique  considered  was  a  receive-adaptive  only  processor  implemented  as 
a  cascade  of  a  normal  adaptive  doppler  processor  and  a  receive-adaptive  polarization 
processor.  This  method  transmits  either  horizontal  or  vertical  fixed  polarization,  while  its 
receive  polarization  adapts  to  minimize  the  clutter  power.  The  doppler  processing  is 
optimum. 

A  fixed-polarization  doppler  processor  is  included  for  reference  purposes.  This  method 
represents  a  system  in  which  both  the  transmit  and  receive  polarization  orientations  are 
either  horizontal  or  vertical.  Again,  the  doppler  processing  is  optimum. 

Finally,  a  noncoherent  processor  is  implemented,  also  for  reference.  Both  the  transmit 
and  receive  polarization  orientations  are  either  fixed  horizontal  or  fixed  vertical. 

A  key  to  the  curves  on  the  PD  plots  identifying  the  processors  associated  with  them  is 
given  in  Table  1. 


TABLE  1.  Key  to  Curves  on  PD  Plots 


Curve 

Transmit 

polarization 

Receive 

polarization 

Doppler 

filterine 

Characteristics 

1 

Pulse-to-pulse 

adaptive 

Pulse-to-pulse 

adaptive 

Adaptive, 

optimum 

Optimum  Matrix 
Method 

2 

Fixed 

Adaptive 

Adaptive, 

optimum 

Receive 

polarization  adapts 
to  minimize 
clutter  power 

3 

Fixed 

Fixed 

Adaptive, 

optimum 

No  polarization 
adaptivity 

4 

Fixed 

Fixed 

None 

Noncoherent 
pulse  summation 

5 

Pulse-to  pulse 

adaptive, 

constrained 

Pulse-to-pulse 

adaptive 

Adaptive, 

optimum 

Constrained  OMM 
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The  processor  evaluation  software  produced  probability  of  detection  (PD)  curves 
using  recorded  target/clutter  data.  PD  plots  were  also  created  for  simulated 
targets  to  predict  radar  performance  under  conditions  for  which  no  recorded  data  was 
available.  In  cases  that  used  simulated  target  or  clutter  inputs,  the  program  was  given 
the  appropriate  parameters.  For  simulated  targets  the  following  parameters  were 
specified: 

a.  Angle  -  This  is  the  primary  orientation  of  the  target  scatterer.  Possible  values 
sire  horizontal,  45°  Unear,  and  vertical. 

b.  Angle  Spread  -  This  is  the  standard  deviation  of  the  angle  of  the  target  with 
respect  to  its  mean  single,  in  radisms. 

c.  Doppler  Spread  -  This  is  the  standard  deviation  of  the  target  doppler 
frequency,  as  a  decimsil  fraction  of  the  pulse  repetition  frequency  (PRF).  The  softwsu-e 
requires  all  targets  to  be  at  zero  mean  doppler  frequency. 

The  psirameters  that  were  selected  in  the  clutter  simulation  are: 

a.  Angle  -  Tnis  is  the  primary  orientation  of  the  clutter  polarization.  Possible 
vsilues  are  horizontsd,  45°  linear,  and  verticsd. 

b.  Angle  Spread  -  This  is  the  stsindard  deviation  of  the  clutter  polarization  single 
with  respect  to  its  mean  angle,  in  radisms. 

c.  Doppler  Mean  -  This  is  the  mesin  clutter  doppler  frequency  as  a  decimal 
fraction  of  the  PRF 

d.  Doppler  Spread  -  This  is  the  standard  deviation  of  the  clutter  doppler 
frequency,  as  a  decimsd  fraction  of  the  PRF. 

These  parsimeters  were  chosen  for  esich  example  to  correspond  to  a  resisonable 
model  for  sin  aircraft  target,  ground  sind  rsun  clutter,  or  chaff,  as  needed  for  the 
simulation  input.  The  output  of  these  test  cases  are  displayed  in  plots  of  PD  versus 
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SNR  for  three  clutter-to-noise  ratio*  (CNR).  The  CNR*  were  -20dB,  -flOdB,  and 
+30dB.  The  complete  set  of  results  is  available  for  several  sets  of  recorded  and 
simulated  data.  The  data  plots  presented  here  represent  the  case  of  recorded 
target /recorded  ground  clutter  with  the  following  target  parameters: 

Angle  =  horizontal 

Angle  Spread  =  0.2 

Doppler  freq  =  0  (normalized) 

Doppler  Spread  =  0.02 

Similar  parameters  apply  for  the  recorded  ground  clutter.  Figures  8  to  13  show 
the  results  of  the  processor  evaluation  for  the  three  clutter  to  noise  ratios  considered. 
Figure  14  is  a  plot  of  the  compressed  target  data  used  in  this  example  when 
transmitting  and  receiving  horizontal  polarization.  The  target  is  at  36  nautical  miles 
in  range  and  is  25  dB  above  noise  and  sidelobe  components.  Figure  15  is  a  plot  of  the 
pulse  compressed  ground  clutter  data  used  in  the  example.  In  this  case,  the  clutter  was 
approximately  20  dB  stronger  them  the  target.  Figure  16  is  a  plot  of  the  puke 
compressed  target  when  the  transmitted  waveform  was  horizontal  and  the  received 
waveform  was  vertical.  The  signal  was  down  approximately  17  dB  from  the  HH  case 
considered  previously.  Figure  17  is  a  plot  of  the  pulse  compressed  clutter  data 
resulting  from  horizontal  transmit  and  vertical  receive  polarizations.  The  clutter  in 
this  case  is  approximately  15  dB  weaker  than  the  HH  case.  Figures  18  and  19  are  the 
target  and  clutter  spectra  for  the  cross  polarized,  horizontal  transmit  and  vertical 
receive  combination. 

The  OMM  and  COMM  processors  showed  superior  detection  results  over  the 
other  processors  for  this  data  set,  as  in  all  the  cases  evaluated. 
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Figure  8.  Horizontal  Transmit  Polarization,  CNR  —  — 20dB 
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Figure  9.  Vertical  Transmit  Polarization,  CNR  =  -20dB 
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Figure  12.  Horizontal  Transmit  Polarization,  CNR  =  30dB 
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Figure  13.  Vertical  Transmit  Polarization,  CNR  =  30dB 
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RANGE  (NMI) 

Figure  16.  Pulse  Compressed  Target  Data  (HV) 
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Figure  17.  Pulse  Compressed  Clutter  Data  (HV) 


4.0  EXPERIMENTAL  RESULTS 


In  this  section,  the  results  of  recent  RADC  polarization  processor  evaluation 
experiments  are  presented.  Actual  flight  test  data  provided  improved  confidence  in  the 
evaluation  compared  to  simulated  data,  and  resulted  in  a  library  of  fundamental 
measurements  of  the  doppler  and  polarization  properties  of  targets  and  chaff. 

The  data  collection,  reduction,  and  processor  evaluation  were  performed  as  previously 
described,  but  results  were  plotted  in  several  different  forms.  Included  here  are  raw  data, 
pulse  compressed  data,  and  spectral  analyses  of  HH,  HV,  VH,  and  VV  components.  In  the 
previous  section,  the  target  and  clutter  data  were  combined  in  the  pre-processing 
algorithm  for  subsequent  evaluation.  In  this  section,  we  present  a  three  dimensional 
range-doppler  plot  of  target  and  chaff  from  the  same  data  recording. 

The  chaff  drop  took  place  in  the  vicinity  of  Griffiss  AFB,  New  York,  in  May  1988, 
using  chaff  cut  specifically  to  enhance  reflections  at  S-Band.  The  target  was  a  medium  size 
aircraft  which  flew  above  the  chaff  cloud.  The  digitized  receiver  outputs  were  recorded 
during  the  chaff  drop,  and  processed  off-line  after  the  flight.  This  experimented  mode 
provides  two  principal  advantages.  First,  the  same  series  of  flight  test  data  was  replayed 
for  processing  by  each  algorithm,  to  insure  a  valid  comparison.  Second,  processing 
algorithms  to  be  evaluated  could  be  arbitrarily  complex  since  they  were  relieved  of  the 
requirement  to  operate  in  real  time.  Hardware  pulse  compression  was  available  for  rapid 
location  of  targets  during  the  test.  This  was  invaluable  for  locating  targets  in  the  mass  of 
recorded  data  during  high  speed  playback.  The  algorithms  employed  software  pulse 
compression  with  floating  point  Fast  Fourier  Transform  (FFT)  computations  to  achieve 
maximum  dynamic  range,  once  the  critical  data  segments  had  been  located. 

Figure  20  is  a  plot  of  uncompressed  data  showing  the  co-polanzation  and 
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cross-polarization  components  from  successive  interleaved  transmit  pulses.  The  transmit 
sequence  allowed  the  four  scattering  matrix  components  to  be  observed  every  two  pulses. 
The  first  pulse  transmitted  was  vertically  polarized,  with  both  vertical  and  horizontal 
receive  data  recorded.  The  next  pulse  transmitted  was  horizontally  polarized,  with  both 
vertical  and  horizontal  returns  recorded  simultaneously.  The  plot  represents  a  range  extent 
of  40  nautical  miles,  4  miles  per  division,  and  a  linear  amplitude  scale.  The  target  can  be 
seen  in  the  VV  plot  at  approximately  12  nautical  miles,  about  4  miles  before  the  edge  of 
the  chaff  cloud.  The  target  also  appears  in  the  remaining  three  plots,  but  at  a  reduced 
amplitude  relative  to  the  background  level.  The  chaff  return  is  maximum  in  the  HH 
channel,  but  has  substantial  components  on  all  plots,  presumably  due  to  atmospheric 
turbulence.  The  signal  components  of  VH  and  HV  are  nearly  equal.  Thermal  noise  is 
independent  between  these  channels,  however,  and  is  slightly  greater  in  the  vertical  receive 
channel. 

Figure  21  shows  data  taken  approximately  three  seconds  after  the  data  shown  in 
Figure  20.  The  fine  structure  of  the  plots  are  somewhat  different,  but  the  signal 
components  have  nearly  the  same  properties  sis  before.  The  chsiff  structure  changes 
continuously  due  to  varying  dipole  orientation  as  affected  by  the  their  fall  rate  and  wind 
velocity.  Figures  22  and  23  are  plots  of  the  compressed  pulse  showing  both  the  target  and 
the  chaff  for  both  of  the  above  csises.  The  experimental  work  described  in  this  paper  was 
performed  in  two  stages.  The  first  effort,  reported  in  reference  [9]  included  eight 
evaluations  using  various  combinations  of  simulated  and  recorded  targets  and  clutter  data 
to  illustrate  detection  characteristics  of  the  processors  operating  both  in  thermal  noise  and 
against  clutter.  Recorded  ground  clutter,  ram  clutter,  and  aircraft  target  data  were  used. 
Simulated  data  used  in  stage  one  included  aircraft  targets,  ground  clutter  and  chaff. 
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VH 

HV 

HH 

Figure  20.  Uncompressed  Target  and  Clutter  Data 
I  Channel  versus  Range 

The  second  stage  of  the  experiment,  conducted  in-house  included  an  upgrade  of  the 
S-Band  radar  sensor  and  adaptive  calibration  of  the  receiver  channels  to  obtaiin  equal  HV 
and  VH  returns.  The  data  analyzed  in  this  stage  was  from  an  actual  chaff  drop  with  am 
aircraft  test  target. 

In  the  first  stage,  the  performance  of  the  candidate  processors  waa  evaduated  for 
several  cases  including  all  combinations  of  recorded  and  simulated  clutter  and  target  data. 
The  processors  that  were  evaluated  include  a  fixed-polarization  noncoherent  integrator,  a 
fixed-polanzation  receive-adaptive  doppler  filter  bank,  the  Optimum  Matrix  Method 
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(OMM),  and  a  Constrained  Optimum  Matrix  Method  (COMM)  processor.  The  algorithms 
were  compared  by  examining  PD  curves  versus  SNR  and  measuring  the  relative 
performance  of  the  processors  expressed  in  dB. 

The  fixed  polarization  noncoherent  integrator  was  used  as  the  basis  for  comparison 
with  the  other  processors.  Probability  of  detection  values  up  to  0.7  were  obtained  for  the 
CNR  =  -20  dB  (thermal  noise  only)  cases  in  which  the  transmitted  polarization  coincided 
with  the  orientation  of  the  target.  PD  values  were  dramatically  lower  for  the  same 
processor  implemented  with  transmitted  polarization  orthogonal  to  the  target.  Cases  using 
recorded  target  and  recorded  clutter  indicated  that  a  12  dB  higher  signal  level  is  required 


to  achieve  the  same  PD  when  the  CNR  is  increased  from  -20  dB  to  +10  dB.  Probability  of 
detection  approached  zero  for  all  cases  where  CNR  =  +30  dB. 

The  fixed-polarization  doppler  processor  models  a  radar  in  which  transmit  and  receive 
channels  have  the  same  polarization,  and  uses  optimum  doppler  processing.  This  processor 
performed  better  than  the  noncoherent  integrator.  Probability  of  detection  was  improved 
over  the  noncoherent  integrator  for  the  CNR  =  -20  dB  (thermal  noise  only)  cases  in  which 
the  transmitted  polarization  coincided  with  the  orientation  of  the  target.  PD  values  were 
dramatically  lower  for  the  processor  implemented  with  transmitted  polarization  orthogonal 
to  the  target.  Cases  using  recorded  target  and  clutter  indicated  that  an  11  dB  higher  signed 
level  is  required  to  achieve  the  same  PD,  when  the  CNR  is  increased  from  -20  dB  to  +10 
dB.  Probability  of  detection  approached  zero  for  all  cases  where  the  interference  was  strong 
(CNR  =  +30  dB). 

The  fixed-polarization  processor  represents  the  conventional  doppler  radar  which  uses 
the  same  linear  polarization  for  both  transmit  and  receive  antennas.  The  mam  cause  of 
poor  performance  for  this  processor  is  its  inability  to  illuminate  targets  whose  primary 
scattering  orientation  is  orthogonal  to  the  radar  receiving  polarization. 

The  receive-adaptive  processor  exhibited  the  worst  detection  capability  of  all  the 
processors.  This  implementation  models  a  radar  system  that  uses  a  fixed  polarization 
transmit  pulse,  but  adapts  the  receive  polarization  so  that  it  is  orthogonal  to  the  clutter 
return.  Doppler  filtering  is  optimum. 

To  achieve  equivalent  detection  performance,  the  receive— adaptive  processor  required 
an  average  of  17  dB  higher  signal  level  than  the  noncoherent  integrator  for  a  CNR  of  -20 
dB,  when  the  transmitted  polarization  coincided  with  the  orientation  of  the  target. 
Detection  was  impossible  when  the  transmit  polarization  was  orthogonad  to  a  target  with  a 
small  angular  spread. 
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Figure  22  Pulse  Compressed  Target  and  Clutter  Data 
Magnitude  versus  Range 

Probability  of  detection  performance  for  the  receive-adaptive  processor  (CNR  =  + 10 
dB)  was  from  2  to  9  dB  poorer  than  for  the  noncoherent  integrator.  Probability  of 
detection  approached  zero  for  all  CNR  =  4-30  dB  cases. 

The  receive— adaptive  processor  represents  a  partially  adaptive  doppler  radar  which 
uses  the  polarization  discriminant  to  achieve  a  receive  polarization  that  nulls  clutter 
returns.  The  experiments  revealed  a  deficiency  in  this  algorithm  in  that  it  sometimes 
severely  attenuates  the  target  along  with  the  clutter.  The  effect  was  observed  in  all  cases  in 
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Figure  23.  Pulse  Compressed  Target  and  Clutter  Data 
3  second  delay  from  Fig  22 

which  the  target  and  clutter  had  the  same  angular  orientation.  As  with  the  other 
fixed-polarization  processors,  performance  was  also  degraded  by  its  inability  to  illuminate 
targets  whose  primary  scattering  orien4  ition  was  orthogonal  to  the  transmit  polarization. 

The  OMM  processor  gave  the  best  detection  performance  of  all  those  evaluated.  This 
processor  models  a  radar  system  that  is  jointly  adaptive  in  both  doppler  and  polarization 
domains.  The  OMM  processor  performed  6  to  12  dB  better  than  the  noncoherent  integrator 
in  the  recorded  target  /recorded  clutter  example  for  a  CNR  of  -20  dB.  PD  values  of  0.95  to 
0.99  were  obtained  at  a  SNR  of  12  dB  in  this  example,  compared  with  values  of  0.05  to  0.55 
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for  the  noncoherent  integrator.  The  OMM  processor  coped  very  well  with  rising  levels  of 
interfering  clutter.  For  a  30  dB  increase  in  CNR,  from  -20  dB  to  +10  dB,  the  OMM 
required  an  average  increase  of  only  9  dB  in  signal  level  to  preserve  detection.  Detection 
performance  for  the  OMM  processor  was  degraded  only  approximately  5  dB  for  CNR  = 
+30  dB  compared  to  that  obtained  for  CNR  =  +10  dB  using  recorded  targets  and  clutter. 

Overall,  the  Optimum  Matrix  Method  processor  demonstrated  far  better  probability 
of  detection  than  the  other  processors,  especially  in  severe  clutter,  where  the  non-optimal 
processors  gave  little  or  no  detection  capability. 

The  Constrained  Optimum  Matrix  Method  (COMM)  processor  demonstrated  the 
second  best  detection  performance  of  the  processors  tested.  This  processing  algorithm  is 
nearly  identical  to  the  OMM  with  the  exception  that  the  transmit  polarization  orientation 
is  selected  from  one  of  16  discrete  angles,  equally  spaced  from  0  to  360  degrees.  This 
represents  a  simplified  implementation  of  the  OMM  algorithm  that  would  require  less 
transmitting  hardware  complexity. 

The  COMM  processor  detection  performance  was  an  average  of  6  dB  worse  than  the 
OMM  implementation  in  the  recorded  target/recorded  clutter  examples  over  all  three  CNR 
levels.  This  processor  shows  a  6  dB  improvement  over  the  noncoherent  integrator. 

The  stage  two  effort  featured  the  availability  of  actual  chaff  and  target  data  for 
further  processor  evaluation.  The  range-doppler  properties  of  the  measured  data  are 
presented  in  Figure  24.  The  doppler  processing  method  consisted  of  applying  40  dB 
Chebychev  weights  to  multiple  pulse  data  from  each  range  index,  then  performing  an  FFT 
to  obtain  the  doppler  spectrum.  The  chaff  appears  highly  spread  in  doppler,  presumably 
due  to  atmospheric  turbulence.  At  S-band,  the  doppler  frequency  is  approximately  12  Hz 
per  knot,  and  doppler  ambiguities  occur  at  multiples  of  the  pulse  repetition  frequency.  The 
target  doppler  is  ambiguous  on  this  plot,  and  has  wrapped  around  several  times. 
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The  results  of  the  chaff  drop  and  target  aircraft  experiments  carried  out  in  the  second 
stage  of  this  effort  show  some  differences  in  processor  performance  compared  to  the  stage 
one  results.  See  Figures  25  through  30.  The  fixed  polarization,  noncoherent  integrator,  haul 
maximum  performance  when  the  transmitted  polarization  coincided  with  the  target 
orientation.  For  the  case  where  CNR  =  -20  dB,  PD  values  up  to  0.5  were  observed.  For 
cross  polarized  returns,  PD  decreased  to  0.15.  Increasing  the  CNR  to  +10  dB  decreased  the 
PD  to  0.05  for  co— polar  and  0.005  for  the  cross-polar  case.  Probability  of  detection 
approaches  zero  for  both  orientations  when  CNR  =  30  dB.  These  results  show  that  the 
noncoherent  integrator  performed  better  in  stage  one  (recorded  target  and  ground  clutter) 
than  m  stage  two  (recorded  aircraft  and  chaff).  The  measured  performance  difference  was 
1.5  dB  for  the  co-polar  return  and  0.5  dB  for  the  cross-polar  return. 

As  in  the  first  stage,  the  fixed  polarized  doppler  processor  performed  better  than  the 
noncoherent  integrator.  In  the  case  of  CNR  =  -20  dB,  the  observed  PD  was  0.55  for  the 
co-polar  return,  and  0.30  for  the  cross-polar  return.  For  the  case  of  CNR  =  10  dB,  the  PD 
became  0.25  for  the  co-polar  return  and  0.05  for  the  cross-polar  return.  PD  dropped  to 
zero  for  CNR  =  +30  dB.  These  results  showed  somewhat  higher  performance  than  the 
stage  one  results,  i.e.,  6  dB  for  the  co-polar  return  and  4  dB  for  the  cross— polar  return, 
with  CNR  =  -20  dB. 

The  fixed  polarization  processor  achieved  a  PD  of  0.65  for  co— polar  returns  and  0.35 
for  the  cross-polar  returns  for  CNR  =  -20  dB.  For  CNR  =  +10  dB,  PD  was  reduced  to 
0.25  for  co-polar  returns,  and  0.05  for  cross-polar  returns,  and  at  CNR  =  +30  dB,  PD 
approached  zero.  As  compared  to  stage  one,  these  results  show  a  decrease  of  4  dB  for 
copolar  returns  and  an  increase  of  3  dB  for  the  cross-polar  returns  when  CNR  =  -20  dB. 
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Figure  24.  Range  Doppler  3D  Plot  of  Aircraft  and  Chaff 

The  receive-adaptive  processor  showed  mixed  results.  While  performance  was  poorer 
than  the  other  processors  with  CNR  =  -20  dB,  it  was  better  than  the  noncoherent 
integrator  and  fixed  polarization  processor  when  CNR  =  +10  dB.  The  PD  for  CNR  =  -20 
dB  was  0.25  for  co-polar  returns  and  0.15  for  cross-polar  returns.  With  CNR  =  +10  dB, 
PD  became  0.15  for  co-polar  returns,  and  0  075  for  cross-polar  returns,  and  for  CNR  = 
+30  dB,  PD  approached  zero.  These  results  represent  an  increase  of  more  than  7  dB  for 
both  co-polar  and  cross-polar  operation  over  stage  one  results,  with  CNR  =  -20  dB. 
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The  OMM  processor  produced  the  best  results,  as  observed  in  stage  one.  For  CNR  = 
-20  dB,  a  PD  of  0.8  was  obtained  for  co-polar  returns,  and  0.7  for  cross-polar  returns. 
With  CNR  =  +10  dB,  the  PD  became  0.55  for  both  co-polar  and  cross-polar  returns.  In 
the  case  of  CNR  =  +30  dB,  PD  dropped  to  0.001  for  both  co-polar  and  cross-polar 
returns.  These  results  were  5  dB  poorer  than  those  observed  in  stage  one,  most  likely  due 
to  the  turbulent  and  noiselike  nature  of  the  chaff  returns  as  compared  with  ground  clutter 
The  COMM  processor  was  second  in  performance  to  the  OMM  processor,  as  in  stage  one, 
achieving  a  PD  of  0.45  for  both  co-polar  and  cross-polar  returns  with  CNR  =  +10  dB. 

Refer  to  Table  1  for  a  complete  key  to  PD  curves  in  Fig  2o  -  Fig  30. 

In  summary,  the  processor  identification  is: 

1  -  Optimum  Matrix  Method 

2  -  Receive-adaptive  polarization 

3  -  Fixed  polarization,  optimum  doppler 

4  -  Fixed  polarization,  noncoherent  pulse  summation 

5  -  Constrained  Optimum  Matrix  Method 
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Figure  27.  Horizontal  Transmit  Polarization,  CNR  =  10  dB 
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Figure  28.  Vertical  Transmit  Polarization,  CNR  =  10  dB 
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5.0  CONCLUSIONS 


The  foregoing  experimental  results  have  demonstrated  that  adaptive  coherent 
polarimetric  radar  techniques  provide  enhanced  target  detection  in  the  presence  of  clutter 
and  chaff.  An  analytic  approach  was  developed  for  dual  channel  optimum  processing  of 
range-doppler  spread  targets  and  clutter.  Matrix  methods  were  utilized  to  develop  and 
implement  an  optimal  polarization  diverse  transmit  waveform  and  optimum  dual  channel 
receiver,  which  maximize  the  signal-to-interference  ratio  and  enhance  target  detection. 
This  doubly  adaptive  approach,  referred  to  as  the  Optimum  Matrix  Method  (OMM),  was 
shown  experimentally  to  be  superior  to  conventional  radar  processing  techniques.  A 
Constrained  OMM,  limited  to  one  of  sixteen  fixed  transmit  polarizations  was  also 
evaluated  and  nearly  approaches  the  performance  of  the  fully  adaptive  system  Three 
conventional,  fixed  polarization  radar  systems  were  implemented  and  compared  to  the 
polarization  diverse  processors.  The  average  performance  improvement  which  could  be 
attributed  to  adaptive  polarization  diversity  was  6-10  dB  for  the  Optimum  Matrix 
Method  and  4-6  dB  for  the  Constrained  OMM. 
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Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (CiI)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  CSI  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability / maintainability  and  compatibility. 
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